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We develop a bosonization scheme for the two-dimensional electron gas in the presence of an 
uniform magnetic field perpendicular to the two-dimensional system when the filling factor is one 
[y = 1). We show that the elementary neutral excitations of this system, known as magnetic 
excitons, can be treated approximately as bosons and we apply the method to the interacting 
system. We show that the Hamiltonian of the fermionic system is mapped into an interacting 
bosonic Hamiltonian, where the dispersion relation of the bosons agrees with previous calculations 
of Kallin and Halperin. The interaction term accounts for the formation of bound states of two- 
bosons. We discuss a possible relation between these excitations and the skyrmion-antiskyrmion 
pair, in analogy with the ferromagnetic Heisenberg model. Finally, we analyze the semiclassical 
limit of the interacting bosonic Hamiltonian and show that the results are in agreement with those 
derived from the model of Sondhi et al. for the quantum Hall skyrmion. 

PACS numbers: 71.70.Di, 73.43.Lp, 73.43.Cd 



I. INTRODUCTION 

Bosonization of fermionic systems is a nonperturbative 
method which has become a very useful tool for treating 
strongly correlated systems in one-dimension. The basic 
idea of this approach consists of describing the neutral 
elementary excitations of the system in a bosonic lan- 
guage, which allows us to map the sometimes intractable 
fermionic system into a more friendly bosonic model. 
A very detailed description of the so-called construc- 
tive one-dimensional bosonization method, its relations 
with the field-theoretical approach and some applications 
can be found ini£ and the references therein. Some ef- 
forts have also been made to extend this formalism to 
higher dimensions. The first attempt was carried out 
by Luthefl£ and extended by Haldanei^ A bosonization 
method for a Fermi liquid in any number of dimensions 
was constructed by Castro Neto and Fradkm 5 and also 
by Houghton and Marston&i 

The quantum Hall effect is one of the most interest- 
ing phenomenon observed in the two-dimensional elec- 
tron systems (for a review see Refs4*2ii2iii). While the 
integral quantum Hall effect can be understood in terms 
of a noninteracting electron model, correlation effects due 
to the Coulomb interaction between the electrons are im- 
portant to understand the fractional quantum Hall effect. 
An exception to the above scenario is the quantum Hall 
system at filling factor one (y = 1), where the electron- 
electron interaction also plays an important role. 

A bosonization approach for the two-dimensional elec- 
tron gas subject to an external perpendicular magnetic 
field was developed by Westfahl Jr. et alH In this case, 
it was shown that the elementary neutral excitations of 
the system, known as magnetic excitons^, can be de- 
scribed in a bosonic language and the Hamiltonian of 



the interacting two-dimensional electron (2DEG) gas was 
mapped into a quadratic bosonic Hamiltonian. However, 
this method can be applied only in the limit of small ex- 
ternal field when a large (integer) number of Landau lev- 
els are completely filled. A different bosonization scheme 
for the collective dynamics of a spinless 2DEG in the low- 
est Landau level was developed by Conti and Vignale. 15 

As pointed out above, the two-dimensional electron gas 
at v = 1 is a strongly correlated electron system. It is 
well known that the ground state of this system is a spin- 
polarized state in which all electrons completely fill the 
lowest Landau level with spin up polarization (quantum 
Hall ferromagnet)i&2ii2iii The elementary neutral exci- 
tations are also described as magnetic excitons^ which, 
in the long wavelength limit, can be seen as spin wave 
excitations of the the quantum Hall ferromagnet. More- 
over, the low lying charged excitation is described by a 
charged spin texture, called quantum Hall skyrmionJ£i 
This non trivial excitation can be viewed as a configu- 
ration in which the spin field at the center points down 
and then it rotates smoothly as we move radially out- 
wards from the center until all the spins point up as in 
the ground state. These charged spin textures are topo- 
logically stable objects with size (the number of reversed 
spins) fixed by the competition between the Coulomb and 
Zeeman interactions. 

Since the quantum Hall system at v = 1 is a quite 
interesting strongly correlated electron system and the 
bosonization was successful in describing the integral 
quantum Hall system at v 3> 1, we would like to extend 
the methodology developed by Westfhal Jr et alml to the 
case when the system is under the effect of a high exter- 
nal magnetic field, in particular, when the filling factor 
is one. This is precisely the aim of this paper. 

Following the ideas of Ref. and0, we start with a 
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Landau level description of this system and then we in- 
troduce a nonperturbative bosonization approach for it. 
We follow Tomonaga's ideas for the one-dimensional elec- 
tron gas in order to show that the neutral excitations of 
2DEG, the electron- hole pairs called magnetic excitons, 
can be treated approximately as bosons. 

We will show that the Hamiltonian of the interacting 
two-dimensional electron gas at v — 1 can be mapped 
into an interacting bosonic Hamiltonian, where the single 
particle energy of the bosons is equal to the energy of the 
magnetic excitons derived by Kallin and Halperin. 13 The 
interaction between the bosons gives rise to the formation 
of two-boson bound-states. In analogy with the isotropic 
Heisenberg model, these bound states can be related to 
the skyrmion-antiskyrmion pair, which is also a neutral 
excitation of the system. 

The paper is organized as follows. In Sec. [HI we will 
present the bosonization scheme for the 2DEG at v = 1, 
i.e, the bosonic operators will be defined. The bosonic 
representation of the density and spin operators will be 
derived and we will show the relation between the low- 
est Landau Level (LLL) projection formalism and the 
bosonization method; the reorganization of the Hilbert 
space will be also discussed. In Sec. IIIII we will ap- 
ply the method to study the interacting two-dimensional 
electron gas at v = 1. Finally, in Sec. IIVI we will an- 
alyze the semiclassical limit of the interacting bosonic 
Hamiltonian derived in the previous section following the 
procedure of Moon et tdmi 



II. THE BOSONIZATION METHOD FOR THE 
2DEG AT v = 1 



A. bosonic operator definition 



Let us consider iV spinless noninteracting electrons 
moving in the x — y plane (two-dimensional electron gas) 
in an external field B = Bqz. The system is described 
by the Hamiltonian 



n = — I d'r *T( r ) (-iftV + -A(r) ) *(r), (1) 



where m* is the effective mass of the electron in the 
host semiconductor (see Appendix lA^I and ^(r) is the 
fcrmionic field operator. In the symmetric gauge, the 
vector potential A(r) = — (r x B)/2 and therefore the 
fermionic field operators can be written in a Landau level 



basis as ( Appendix lB|l 

* f ( r ) = X)(»»»|r}cJ, m 

n,m 

- V - e- |r|2/4/2 G* (ir/Dc 1 

*(r) = J^(r|n m)c nm 

n,m 

- V 1 e~^ 2/Al2 G . (ir/Dc 
t Z V 2ttZ 2 



(2) 



where r = x + iy, I = y^hc/ (eB) is the magnetic length 
and the function G m +n,n(x) is defined in Appendix IU1 
The fermionic operator c\ lm creates an electron in the 
Landau level n with guiding center m and obeys canoni- 
cal anti-commutation relations 



{ c nmi c n'm'} — { c nrm c n'm'} 



(3) 



with n = 0, 1, 2, . . . and m = 0, 1, . . . , N^ — l. Here, N# = 
AriB is the degeneracy of each Landau level, with rig — 
1/(2ttZ 2 ) and A is the area of the system. Substituting 
Eqs. J5J in Eq. QJ, we find that the Hamiltonian of the 
system is diagonal in the Landau level basis, 

1, 
2' 



U Q = ^2hw c (n + ^)c\ m c nm , 



(4) 



where w c = eB/(m*c) is the cyclotron frequency. Defin- 
ing the filling factor v = N/N^ as the number of filled 
Landau levels, for v integer, the ground state of the two- 
dimensional electron gas (2DEG) is obtained by com- 
pletely filling the v lowest Landau levels, 



\° s ) = n n c «™i°) 



(5) 



n— m— 



where |0) is the vacuum state. 

Now, if we consider the electronic spin and restrict the 
Hilbert space to the lowest Landau level (n — 0) only, 
the fermionic field operators J2J) become 



(6) 



a i 



where c m a creates a spin a electron in the lowest Landau 
level with guiding center m. These creation and annihila- 
tion fermionic operators also obey the anticommutation 
relations, 



{ c mcn c m' a'} ~ { c m a, Cm 1 a'} — 0, 



(7) 
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FIG. 1: Schematic representation of (a) the 2DEG ground 
state at v = 1 (quantum Hall ferromagnet) and (b) the ele- 
mentary neutral excitation (the electron-hole pair), g is the 
Zeeman energy. 



with a =T or J.. 

In addition to the kinetic energy term Ho [Eq. JIJ], 
we should also include a Zeeman energy term Hz in the 
total Hamiltonian of the system, 



(8) 



where g* > is the effective electron g-factor (see Ap- 
pendix [3J and hb is the Bohr magneton. Substituting 
Eqs. © in the expressions Q and JSJl, the total Hamil- 
tonian of the 2DEG, H = H + H z , yields 



-g*H B B ^ 



°~ C ma c mcr- (9) 



We can see that H is also diagonal in the Landau level ba- 
sis and the kinetic energy term is simply a constant. The 
one particle energy eigenvalues are — g*/xs/2 and g* //b/2 
whereas the degeneracy of each energy eigenstates is N^. 

For v = 1, the ground state of the 2DEG, \FM), is 
obtained by completely filling the spin up lowest Landau 
level (the quantum Hall ferromagnet) 



\fm)= n 



10), 



(10) 



m=0 



as illustrated in Fig. QJa) . In this case, each guiding cen- 
ter is occupied by only one electron with spin up. \FM) 
is an eigenvector of the operator S z (the z component of 
the total spin) whose eigenvalue is equal to N$/2. 

The neutral elementary excitations of the system are 
electron-hole pairs or spin flips as one spin up electron is 
annihilated and one spin down electron is created in the 
quantum Hall ferromagnet [Fig. 01(b)]. These excited 
states \^f) can be constructed by acting with the spin 
operator S~ — S x ~ iS y on the ground state \FM), 

I*) oc S~\FM). 



In the bosonization approach for the one-dimensional 
electron gas, the annihilation and creation bosonic oper- 
ators are derived from the electron density operator p(k) 
as the electron-hole pair excitations can be obtained by 
acting with p(k) on the ground state of the systemic. In 
order to define the bosonic operators, the commutation 
relation between the electron density operators p(k) with 
different momenta is analyzed. We will follow the same 
procedure, but here we will analyze the spin density oper- 
ators S+(r) = S x (r)+iS y (r) and S~(r) = S x (r)-iS y (r) 
in order to define the bosonic operators for the 2DEG at 
v = 1. More precisely, we are interested in the Fourier 
transform of these spin operators. 

Before doing that, we need to say some words about 
the density operator of spin a electrons, which is defined 
as 

M r ) = *t(r)¥ a (r). 

With the aid of the expressions © and the definition 
of the function G m . m '(x) (Appendix Q , it is possible to 
show that 



Per(0[) 



r e 



= E /rf 2 re- lq - r (m|r)(r|m') C t_c m ^ 

m,m' 

m,m' 

— e I ^ / ^ ^ G m . m > (Iq^c^ fjCyn' (j, (11) 



where q — q x + iq y and the operator z is defined by the 
Eq. I|C1(I . Observe that the action of p a (q) on \FM) 
does not create any electron-hole pair excitations. 
The spin density operator is defined as (ft = 1) 

where the components of the vector a are the Pauli ma- 
trices. However, we will define the spin density operators 
S l+ (r) and ^"(r) only by the nonzero matrix element. In 
terms of the fermionic field operators we have 

S+(r) = *}(r)*i(r), 
S~(t) = *J(r)* T (r). 

In analogy with equation it is easy to show that the 
Fourier transform of these spin operators is given by 



S+(q) 
S-(q) 



-\lq\ 2 /2 



-\lq\ 2 /2 



f~] C*m,m'il<l)c m ^C' 



(12) 



2J G m . m '{lq)c s m[ c m ' t- (13) 
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After some algebra, it is possible to show that the 
commutation relation between the operators 5 + (q) and 
S~(q) is proportional to the Fourier transform of the 
density operators p-\ (r) and pi (r) (see Appendix [D| , 



J 2 qq"/2 



/5 T (q + q') 



Pi(q + q') 



Now, as the average values of /3j(q) and /3j.(q) in the 
ground state (fTTTft are (/5-f (q + q')) = A^<5 q+ q',o and 
(Pi(q + q')) = 0' respectively, the average value of the 
commutator S~,] is 



([^ + ; ^]>=^ |i9|2/2 ^ 



q+q',o- 



(14) 



The above expression will allow us to define the bosonic 
operators as a function of the fermionic operators c] na 
and c m (j . 

If we define the operators 6 q and 6 q by 



6 q = 



6i 



^ e " 9|2/4 5 a , 
/AU q 



(15) 
(16) 



with n q > and ^2 ^q < A^,. 

The state b\\FM) is a linear combination of electron- 
hole excited states as illustrated in Fig. ^b), where both 
the electron and the hole have spin down. In fact, the 
bosonic operator 6 q is similar to the operator e^^q), 
with n = p = 0, discussed in Ref. [T^- This operator cre- 
ates the neutral excitations known as magnetic excitons 
when it is applied on the ground state of the noninter- 
acting two-dimensional electron gas [Eq.JSJ)]. 

The momentum q is canonically conjugate to the vec- 
tor R = (R§ + Ro)/2. Here, the vectors Rg and R§ 
are, respectively, the position of the guiding centers of 
the electron and the hole excited in the system as de- 
fined in Appendix [Bj Hence R is the center of mass of 
the guiding centers of the excited electron and hole. In 
addition, the momentum q is a good quantum number 
because the total momentum of a two-dimensional sys- 
tem of charged particles in an external magnetic field is 
conserved when the total charge of the system is zero^£ 

To sum up, we can say that the state 6 q |FM) is a neu- 
tral elementary excitation of the 2DEG at v = 1 which 
corresponds to either a spin-flip or a magnetic exciton 
with momentum q. 



and if we approximate the commutation relation between 
the bq and & q by its average value in the ground state 

cm, 



K,bl]^(K,bl]) 



(17) 



we can say that 6 q and 6 q are approximately bosonic 
operators. In analogy with the Tomonaga's model for 
one-dimensional electron gas, we will assume that (|17|l is 
exact. 19,20 This is the main approximation of our method. 
Observe that this approximation is quite similar to the 
one adopted in the random phase approximation. 21 
Therefore, after this point, we will assume that 



1 



&q 



=e -M74 y q 

^ ra.m' 



'(-'9)4, 



(18) 



'(i?)4a<WT» (19) 



are bosonic operators, which obey the canonical commu- 
tation relations 



[blbl] = [&„&,/] =0. 



(20) 



The quantum Hall ferromagnet \FM) is the vacuum 
state for the bosons as the action of the fermionic opera- 
tor c m i i on \FM) is equal to zero. Therefore the bosonic 
Hilbert space is spanned by applying the operator & q on 
the quantum Hall ferromagnet any number of times 



iw)= n ^%io>=n 1 %=i^>, 



(21) 



B. Density operator 

Although the bosonic operators are not directly de- 
rived from the electron density operator as in the one- 
dimensional electron gas, the latter is a very useful oper- 
ator when the Coulomb interaction between the electrons 
of the 2DEG is considered. In this section, we will show 
that it is possible to write down the electron density op- 
erator as a product of the bosonic operators 6 q and fe q . 

The bosonic representation of any operator O — 
Ofcjn (T , Cm^cr) can be obtained applying this operator 
on the eigenstates (j2U, which span the bosonic Hilbert 
space, 



oiw> = o n 



(&t)«. 



\FM) 



(22) 



[0,H^L=]\FM) 



^=0\FM). 



Starting with the expressions of O and 6 q as a function 
of the fermionic operators cS m a and c m ' a , we can cal- 
culate the value of 0\FM) as well as the commutation 
relation [0, 6 q ], which allows us to obtain the value of 
the commutator in Ea. (|22|l . 

Following the above procedure, let us derive the expres- 
sion of the density operator of spin up electrons p-\ (k) as 
a function of the b's. It is quite easy to show that 



[/3 T (k),6 q ] 

where k A q = l 2 z 
of the function G m m 



|/fc| 2 /4g-ikAq/2^t 



k+q' 



(k x q). Using the property ljC12|) 
(x), we can see that the value of 
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Pl(k)\FM) is simply a constant, 
p T (k)|FM) = e -l (fe l 2 / 2 G m , m ,(/fc)4 TCm - T |FM) 



= ^.o^lFM). 
After some algebra, we end up with 



(23) 



Mk)|{n q }> 



V p -|ifc| 2 /4-*kAp/2,t 

e °k+p 

p£{n q } 



( 6 t)«p 



x II 

9£{«q},<?#P 



|FAf> 



+ (5k,o^|{n q }) 

= -e- m2/4 Y,^ kAP/2b U p b P \M) 
p 

In the second step above, we used the fact that 



(24) 



>(bl) 



(blT*b r 



and that b p \FM) = 0. As Eq. is valid for any eigen- 
state of the bosonic Hilbert space, i.e., it is an operator 
identity, we can conclude that 



/5|(k) = (5k,o^0 - e 



-\lk\ 2 /4 



£ e - ik A q /2 6 t +q6q (25) 



is the bosonic representation of the density operator 
PT0O- 

In the same way, the expression of the operator pi (k) 
in terms of the 6's is given by 



Pi(k) 



-|ifc| 2 /4 



>p -MkAq/2. t 5 



(26) 



as pi(k)\FM) = and 

[p i (k),6 q]= e-l*l 2 /4 e +^A q /2 6 t +(j . 

An alternative procedure to obtain Eqs. (|25|l and i|26|) 
consists of looking for an expression in terms of the 6's 
which satisfies the commutation relation [O, & q ]. For in- 
stance, for the electron density operator /3f(k), the com- 
mutator [y5|(k), 6 q ] oc fr q+k and therefore we can conclude 
that the expression of p j (k) in terms of the 6's should be 
a linear combination of the operator & q+k & q - Choosing 
the coefficients properly, we easily find the first term of 
Ea. (|25H . In order to obtain the complete expression, it is 
necessary to add the term related to the action of Pf (k) 
on \FM). In the next sections, we will adopt this proce- 
dure as it is simpler than the one previously discussed. 



Finally, from equations (|25|l and l|26(l , we arrive at the 
bosonic form of the density operator 

p k = p T (k)+ Pi (k) (27) 

= <5k,o^0 + 2ze-' ifc l 2 / 4 J2 sin ( k A q/ 2 )4+qV 
q 

Notice that pk is quadratic in the bosonic operators. 

C. Spin density operators 

In this section, we will derive the bosonic representa- 
tion of the spin operators »S k , S k and S^. We will see 
that the obtained forms for these operators are similar to 
the ones of the formalism suggested by Dyson to study 
spin waves in a ferromagnetic system^ 

Since the z-component of the spin density operator can 
be defined as 

S"(r) = i(*{(r)* T (r)-*I(r)* i (r)) 1 

the expression of the Fourier transform of this operator as 
a function of the bosonic operators follows immediately 
from Eqs. and (J5BJ, 



Si, 



1 



(p T (k)-p;(k)) 



(28) 



' /fc ' 2 / 4 ]Tcos(kAq/2)& k+q V 



In the spite of defining the bosonic operators 6 q and 6 q 
from equations l|15|) and respectively, the latter do 
not correspond to the bosonic representation of and 



5 k as 



In fact, we should also follow the procedure described 
in Sec III Bl in order to calculate the bosonic form of the 
operators S£ and . 

From equations i|13|) and l|19|l it is possible to show 
that [S^j&q] = 0. Therefore, the action of 5 k on the 
eigenstates (|21|l is related to the action of this operator 
on the quantum Hall ferromagnet \FM) only 



S k |{« q }) 







1-r (&£)™ q 

+ [] \^=S^\FM). 



?6{» q } V "1- 

Unlike the results for the electron density opera- 
tor, S^lFM) is not just a constant, but rather 
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proportional to a linear combination of the terms 
G m , m i{lq)(y m ,c m /^\FM). In this case, it seems quite rea- 
sonable to consider Eg. l|16|l and define the bosonic rep- 
resentation of the operator as 

5 k = ^e^^bl (29) 
In the next section, we will show that (|29|l is very well 



defined because it satisfies the commutation relations be- 
tween the spin density operators as well as between the 
spin density and electron density operators. 

On the other hand, S£\FM) = [see Ea.ljP2)l] and 
hence the bosonic representation of the operator S£ is 
completely determined by the commutation relation be- 
tween this spin operator and b]., 



[s+,b{] = ^( e ^*/ 2pT(k + q) _/ fe *,/ 2pi(k + q) ) 



,/ l,A( q -p)- q A^ CA 



(30) 



As the first term of Ea. (|3()(l is proportional to (>k,-q, we can conclude that the operator should present a term 

linear in 6k- Moreover, the second term is a linear combination of products of the form & k+p 6 P , which implies that S£ 

should have a term proportional to the product Wbb. Choosing the coefficients properly, we end up with the bosonic 
form of the operator S£ 



5', 



\lk\ 2 /4 



kA(q-p)-qAp ^ frt 



k+q+p^P^q' 



(31) 



which satisfies the commutation relation (|30|l . 

As pointed out earlier, the representation l|28|l . I|29(l 
and (|31|> is similar to the one previously considered by 
Dyson.— An important point of this formalism is that, 
althougth the Hermiticity requirement S£ = (SZ^) does 
not hold, the usual commutation relations between the 
spin operators are satisfied. A detailed review of this 
formalism can be found in Ref. |2J. As we will show in 
the next section, the representation Ij28|l . (|29|l and (|3 1|) 
derived using the bosonization method also preserves the 
commutation relation between the spin density opera- 
tors. 



D. LLL projection 

In this section, we will show that, using the bosonic 
representation of the operators /5k, S£, >5 k and 5 k , the 
commutation relations between them are in agreement 
with the results derived from the formalism of the Lowest 
Landau Level (LLL) projection. 

The LLL projection is a formulation of the quantum 
mechanics in the restricted subspace formed by the lowest 
Landau level as developed by Girvin and Jacb— (a brief 
review of this formalism is presented in Refs.?' 10 ). An 
important consequence of the projection of the electron 
density and spin density operators on the LLL subspace is 
that the commutation relations between those operators 
are modified, i.e., the projected spin operators do not 



commute with the electron density operator and do not 
follow the canonical commutation relations between spin 
operators either. 

From equation (|27|) . it is quite easy to show that the 
commutation relation between electron density operators 
with distinct momenta is given by 

[A,Pq] = 4/ ik '" 2 sin(qAk/2)e-l" +, «l'' 4 

p v ' 

We can see that [p"k, pq\ is proportional to a linear combi- 
nation of the product 6 k+q+p fe p with coefficients equal to 
sin((k + q) Ap/2). This result indicates that the commu- 
tator should be related to the electron density operator 
p q+ k- In fact, if we compare l|32|) with Eq.J37|), we find 
that 

[pk,p q ] = 2ie' 2k ' q / 2 sin (k A q/2) p q+k , (33) 

which agrees with the result obtained from the LLL pro- 
jection formalism. In the LLL projection approach, it is 
proved that the projected electron density operators with 
different momenta obey an algebra similar to the one of 
the translation operators in a magnetic field. When a 
particle in a magnetic field is translated along the paral- 
lelogram generated by the vectors kZ 2 and ql 2 , the parti- 
cle acquires a phase equal to q A k. As a consequence of 
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that, the commutator [pk,Pq] is not equal to zero, con- 
trary to the behavior of the non-projected operators^. 
In the same way, from expressions (|27l) and H28[). we 



find that the commutator between />k and is also dif- 
ferent from zero, 



[p k , S*] = 2ie 1 ' k ' q / 2 sin(k A q/2)e^ lk+l ^ / 4 ^ cos ( v " ■ ^ ' ) b 



(k + q) A p 



] b 

k+q+p°P 



^k,o<5 q , iVj - 2i/ k ' q / 2 sin (k A q/2) S } 



k+q- 



(34) 



This result implies that, within the LLL subspace, the charge and spin excitations are entangled^. As it will be 
discussed in Sees. IIII CI and IIVI the charged excitation of the interacting two-dimensional electron gas at v — 1 is 
described by a charge spin texture (quantum Hall skyrmion)^ 

Finally, after some algebra, it is possible to show that the commutation relations between the spin operators S^, 

S£ and S~ are given by 



[^k ' ^q] 



e l k q / 2 cos(k A q/2) v ^V0e" l ' fe+ ' 9 ' /4 fo k+q 
= / k -<i/ 2 cos(kAq/2)S k+q 
[St,SZ] = -e i2k ^/ 2 cos(kAq/2) 



(35) 



\lk+lq\ 2 /4 



^cosf (k + q + P)A(p - p/) ) b t 



/J^ e -l»+»«l a /4 6 _ t _ q _ 
= -e pk - t i/ 2 cos(kAq/2)5+ +q 
S~] = N,e-^ 2 ^-K - 2e-\^-^ cos (k A q/2) £ cos f^+f^) & t +q+p6p 



'k+q+p+p'kpV 



+ 2 e -N 2 / 4 -l^l 2 /4 sin(kAq/2)Esin ( (k±^Ap\ & t +q+pfop 



2e /2k - q / 2 cos(kAq/2)S, 



k+q 



/ k - q / 2 sin(kAq/2)p k+q . 



r 



(36) 



(37) 



Again, all the commutation relations (|35|l - (|37|l are in 
agreement with the results calculated in the LLL projec- 
tion formalism. 

It is not surprising that our bosonization approach for 
the 2DEG at v = 1 recovers the results obtained with 
the LLL projection. Remember that all operators con- 
sidered until this moment were expanded in terms of 
the fermionic creation and annihilation operators c' m a 
and c m a with the aid of expressions ©, which are the 
fermionic field operators projected in the LLL. In addi- 
tion, as discussed in details in the appendix IO the func- 
tion G m ^ m i(x) is the matrix element in the lowest Lan- 
dau level basis of the projected operator e _lq r . When 
the Fourier transform of any operator is calculated using 
the LLL projection formalism, it is necessary to consider 



the expression of the projected operator e~ lqr . 

Returning to the previous section, we can also conclude 
that the operator S~ is very well defined by Eq.(J25J| as 
this one preserves the commutation relation between the 
electron density and spin density operators within the 
LLL subspace. 



E. Hilbert space 

In the bosonization approach for the one-dimensional 
electron gas, HaldaneSi proved that the Hilbert space 
spanned by an arbitrary combination of the fermionic 
creation and annihilation operators acting on the vacuum 
state | A = 0)o is identical to the Hilbert space spanned 
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by an arbitrary combination of the bosonic creation op- 
erators acting on the set of all TV-particle ground states 
\N)o, with N E Z, which is the vacuum state for the 
bosonsi 

The above assumption can be elegantly proved by 
calculating the grand canonical partition functions of 
the noninteracting fermionic and bosonic Hamiltonians, 
where the latter is derived from the former using the 
bosonization method for the one-dimensional electron 
gas. Since all terms of the partition function are positive 
quantities, the relation between the two functions allows 
us to compare the degree of degeneracy of the fermionic 
and bosonic Hilbert spaces. For the 2DEG at v — 1, we 
have been considering a system constituted by a fixed 
number of particles N = N ( f >1 therefore we will analyze 
the canonical partition function. 

In Sec III Al we showed that the Hamiltonian of the 
2DEG at v = 1 is given only by the Zeeman term [see 
Eq. (O], which is diagonal in the Landau level basis. The 
energy eigenvalues can be written as E n = ng — 
where < n < is the number of electrons with spin 
down. The degeneracy V F of each energy eigenstate can 
be easily calculated, 



TV 



n j \n j \n\(N,f, — n)! / 
Hence, the fermionic partition function is given by 



Z F 



Tr(e 



E 

71=0 



-n(3g 



(38) 



vath.p = l/(K B T). 

On the other hand, as it will be discussed in Sec lIII Al 
the bosonic Hamiltonian [Eq. i|44[l ] derived from the non- 
interacting fermionic one [Eq. (JjJJ] using the bosonization 
scheme is diagonal in the basis of the eigenstates |{n q }) 
[Eq. J2U]. Therefore, the canonical partition function is 
simply given by 

Z» = Tr( e -^)^ ({nq} | e -^| {nq}) 

{«q} 



N,, 



B e~ nl39 . 



(39) 



Here V B is the number of eigenstates with rt-bosons and 
n = n q- Notice that the eigenstates with n > are 
not included in the above sum as those states correspond 
to a number of electron-hole excitations greater than the 
number of fermions N^. 

The values of V B are determined by the number of 
points in the momentum space. The maximum momen- 
tum value can be estimated if we remember that a boson 
of momentum q, created by the action of the operator 
on the state \FM), can be described as an electron-hole 
pair whose distance between the center of their guiding 
centers is |r| = / 2 |q|. Moreover, as discussed in Appendix 



IbI a particle in the lowest Landau level with guiding cen- 
ter m corresponds to a particle moving in a cyclotron 
orbit with radius equal to the magnetic length I whose 
guiding center is located at a distance l\J1m + 1 from the 
origin of the coordinate system. Therefore the largest dis- 
tance between the electron and the hole in the magnetic 
exciton corresponds to m = and it is roughly equal to 
y/2Nj,l. Since the momentum cutoff is q m ax — y2iVj/Z, 
the number of points in the momentum space is given by 



x A f , 2 2nl 2 N <j> f q . 

1 = — j d 2 q= —rj^- I qdq / d9 = N? 



4tt 2 



4tt 2 



(40) 

where A is the system area. 

From the above analysis, the number of states with n 
bosons is given by 





II 


V B 


= Nl 







{Nl-2)\2\ 



1 n 



N'(N' + l)...(N' + (n-l)), n>l, 



hence the canonical partition function for the bosonic 
Hamiltonian can be written as 

Z B = e PgNj2 {1 + g V B e -nfigy (41) 



Since Z B 3> Zq, we can conclude that the bosonic 
Hilbert space is bigger than the femionic one. Even 
having removed the states with n > from the par- 
tition function l)4ip. we still have unphysical states in 
the bosonic Hilbert space. 

There is only one fermionic and one bosonic subspaces 
of the corresponding Hilbert spaces which are identical. 
Notice that the first two terms of equations l|38|l and 
(|41|l are equal, which implies that the fermionic subspace 
spanned by the quantum Hall ferromagnet, \FM), and 
the states with only one spin down (n = 1) is identical 
to the bosonic Hilbert space spanned by the vacuum and 
the states of one-boson b^FM). 

This overcompleteness of the bosonic Hilbert space can 
be easily understood. From expressions l|38|) and (|41|) , we 
can see that the number of states of two-bosons is roughly 
twice the number of fermionic states with two spin down 
(n = 2). If we consider, for instance, a state of two- 
bosons constituted by two bosons of momenta qi and 
q2, such that \lq\\, \lq%\ < 1, in the fermionic language, it 
can be seen as in Fig. [21a. Notice that this state is equiv- 
alent to the one, which is constituted by two bosons of 
momenta q3 and q4, such that \lqs\, {Iq^l 3> 1 [Fig. [21b]. 
Based on that, we can say that we are "double counting" 
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FIG. 2: Schematic representation of a two-bosons state: (a) 
|2gi|, \lqi\ < 1 and (b) \lq 3 \, \lq^\ > 1. 



III. INTERACTING TWO-DIMENSIONAL 
ELECTRON GAS AT v = 1 



In this section, we will apply the bosonization method 
developed for the 2DEG at v — 1 to study the interacting 
two-dimensional electron gas aiv= 1. We will show that 
the Hamiltonian of this interacting system is mapped into 
an interacting bosonic model. 



A. Noninteracting electron system 

As pointed out in the last section, the Hamiltonian of 
the noninteracting two-dimensional electrons at v = 1, 
restricted to the lowest Landau level subspace, is given 
by the Zeeman term only [Eq. @]. In the Landau level 
basis, it can be written as 



the number of states of two-bosons. Of course, this prob- 
lem becomes worse as we take into account states of n 
bosons (n > 2). Besides that, we can still study the low- 
energy physics of the system as, in this case, we have a 
small number of bosons with momentum \lq\ < 1 present 
in the system. As we will see in the next section, the 
energy of the bosons increases with the momentum [Eq. 
(|51|l ] . It is worth mentioning that this same problem ap- 
pears in a description of a bilayer quantum Hall system at 
total filling factor one (i>t = I)i2£ Here, the sponaneous 
interlayer phase coherent (111) state can be viewed as 
a condensate of interlayer particle- hole pairs (excitons), 
which, in the very dilute regime, can be treated as point- 
like bosons. The corresponding bosonic Hilbert space is 
also overcomplete. 

This problem could be fixed, for instance, introducing 
a constraint which restricts the bosonic Hilbert space to 
the physical states only. However, this is quite a hard 
task. For example, it is not possible to follow the ideas of 
the well-known expansions of the spin operators in terms 
of bosons, such as the Schwinger boson representation. 28 
In this case, the local spin operators are written as a 
function of the local bosons operators a.; and a.; , namely 



S; = a]c 



Sf = (a}ai - a}ai)/2. 



The constraint is easily determined as it is related to the 
fact that the number of bosons on the site i should be 
twice the spin S, i.e., a\ai — aja,i = 2S. The same idea 
can not be applied to our case as the bosonic operators 
and bq are not local. In fact, they involve a linear combi- 
nation of electron- hole excitations where the particles are 
localized in different guiding centers. Until now, we have 
not found a systematic way of introducing a constraint 
in our formalism. 



7~Lo = 7~Lz 



'EE' 



(42) 



where g — g* [i^B > 0. 

In order to find out the bosonic form of the Hamilto- 
nian H42JI . it is necessary to calculate the commutation 
relation between 7io and the bosonic creation operator 



*E E 



cr m.n.n' 



-\lq\ 2 /4 



-<rG n ,n'(k) 



X [ c m cr c ™' CT) c n l c n' t] 



5 ^ 1 e -|^| 2 /4 Gnn , ( ^ )c t , T 
n.n' V "r 



(43) 



Since the above commutator is proportional to 6^, Hq 
should present a term of the form g ^ q fr^feq, which gives 
the same commutation relation as in Eq. 1)43(1 . Moreover, 
the action of Hq on \FM) is equal to a constant, —gN^/2. 
Therefore we can conclude that the bosonic form of the 
Zeeman term is 



1 



(44) 



The same result can be obtained in a more rigorous way 
by explicitly calculating the action of Tio on the eigen- 
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states ijHJ, 



bosonic system, whose dispersion relation is constant. 



Wo|{n q }) = H ([[ y -^=\FM)) 



/riq! 



(45) 



[Wo, n 4M fm > + n -^=n \FM) 

\/n a \ \/n a \ 

q V 4 q V m 

(.9 e ^--^n^i™). 

fcG{n q } q V n q- 



This analysis shows that the Hamiltonian of the nonin- 
teracting two-dimensional electron gas at v — 1, restrict 
to the lowest Landau level is recast in a noninteracting 



B. Interacting electron system 

Now, we will consider an interacting two-dimensional 
electron gas at v = 1, where the particles are restricted to 
the lowest Landau level. The Hamiltonian of the system 



H = Ti-o + TLi 



(46) 



Here, H$ is given by Eq. (|42[) and the interacting term 
can be written as 



Hint = \Y,J Mr ' *t(r)^(r')^(|r-r'|)^(r')* CT (r), 



(47) 



where U(|r|) = e 2 /(e?*) is the Coulomb potential and 
e the dielectric constant of the host semiconductor (see 
Appendix lAl . Substituting Eq.© in Hi ntl it is possible 
to write down the interacting term as a function of the 
density operators of electrons a as 



H mt = \Y,j n*)P«r(k)M-k) 



(48) 



where V(k) is the Fourier transform of the Coulomb po- 
tential in 2D, 



V(k) 



2-rte 2 



ek ' 



and k = |k|. Using the bosonic form of the density oper- 
ators /^(k), we can write down the interacting term as 
a function of the bosonic creation and annihilation oper- 
ators. Substituting Eqs. and (UHl in Eq. l|3g)l. we 
have four distinct terms, Hint = H\+H2+Hz + Ha- The 
first one is a constant related to the positive background, 

Ul = i/ y(k)<5k > ' 
whereas the second and third terms are equal to zero as 

H 2 = -H 3 = -t^E [d'kVWN^.oe-W 2 / 4 
71 p J 

x sin(kAp/2)fo q & q . 

The last term is quartic in the bosonic operators. Rewrit- 
ing H4 in normal-ordering in the operators 6, we end up 



r 



with a quadratic and a quartic terms in the bosonic op- 
erators, namely 



H4 — 



0^2 / d 2 ^(k)e H ' fc|2/2 sin(kAp/2) 
n p,q •* 

x sin(k A q/2)6] c+q 6 q 6^_ k 6p 

= / rf 2 ^(k)e-" fc ' 2 / 2 sin 2 (kAq/2)6 q fo q 

n q J 

+ JL^ j d 2 fcl/(k)e-l ifc l 2 / 2 sin(kAp/2) 
n p,q J 

x sin(k A q/2)6] c+q 6 p _ k 6 q & p 

= jyfzo-^'^odMv^^q 

q 

+ / rf 2 fcnk)e- |ifc|2/2 sin(kAp/2) 

p,q 

x sin(k A q/2)6 k+q 6 p _ k 6 q & p . (49) 



Here Io(x) is the modified Bessel function of the first 
kind^ 

Now, adding the bosonic form of the noninteracting 
Hamiltonian [Eq. I|44(l ] to Hint, we will arrive at the ex- 
pression of the total Hamiltonian of the interacting elec- 
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1.25 I , 1 , 1 , 1 , 1 r 




FIG. 3: Dispersion relation of the bosons [Eq. I51L with 
g = 0], in units of the Coulomb energy e 2 /(d) as a function 
of the momentum q. 

trons as a function of the bosonic operators, 

q 

+ 4 E ^(fc)e-l Jfc l 2 / 2 sin(kAp/2) 

k,p,q 

xsin(kAq/2)6] c+q 6t _ k fe q 6 p; ( 50 ) 
where the dispersion relation of the bosons is given by 

Wq = 9 + Tl{l ( l eHi9 ' 2/4/ o(NI 2 /4)) • (51) 

This curve is plotted in Fig. [3] for the case g — 0. The 
energy of the bosons is given in units of the Coulomb 
energy e 2 /(el). 

The Hamiltonian (|5()|l describes a system of interacting 
two-dimensional bosons. The ground state is the vacuum 
state \FM) [Eq. and its energy is equal to Eq = 

— ^gN,/,. This result implies that the ground state of the 
two-dimensional electron gas at v = 1 does not change 
as we moved from the noninteracting to the interacting 
system. As pointed out earlier, the ground state of the 
interacting two-dimensional electron gas is the quantum 
Hall ferromagnet even in the limit of vanishing Zccman 
energy (g -> 0). 

The elementary neutral excitations are described by 
the bosons 6, whose dispersion relation w q is equal to 
the previous diagrammatic calculations of Kallin and 
Halperini 3 . and the results of Bychkov et alml The long 
wavelength excitations can be considered as the spin wave 
excitations of the quantum Hall ferromagnet, while the 
ones with large momenta correspond to a quasielectron- 
quasihole pair as the particles are very far apart. Remem- 
ber that, as discussed in Sec. Ill Al the distance between 
the guiding centers of the excited electron and hole is 
|r|=Z 2 |q|. 



The above results corroborate our initial association 
[Sec. Ill Aj between the bosons b and the magnetic exci- 
tons described in Refll3l As a matter of fact, within our 
bosonization method, we goes beyond the diagrammatic 
calculation^ as we end up with an interaction between 
the bosons. In the next section, we will study the result- 
ing interacting bosonic model. 

C. Bound states of two-bosons 

In this section, we will study the states of two-bosons. 
More precisely, we will check if the interacting bosonic 
model (|50|) can describe the formation of bound states of 
two-bosons. Our initial motivation is based on previous 
results of the one-dimensional ferromagnetic Heisenberg 
model. 

The problem of two interacting spin waves (magnons) 
in the ferromagnetic Heisenberg model was analyzed by 
Dyson 2 ^, who derived a bound state condition when the 
total momentum of the pair is equal to zero. For an 
arbitrary value of the spin S, it was verified that this 
condition is not fulfilled in two and three dimensions and 
it was concluded that bound states of two spin waves 
do not exist, in contradiction to the results for the one- 
dimensional models After that, this problem was also 
discussed by Wortis^i who, in opposition to Dyson's re- 
sults, proved the existence of bound states of two spin 
waves for any value of the spin S and the dimensionality 
of t he sy stem. A review about this topic is presented in 
Ref. Hi 

On the other hand, Tjon and Wright 32 studied the 
dynamical solitons of the one-dimensional ferromagnetic 
Heisenberg model. The dynamical soliton is a solution 
of the dynamical equations of motion, localized in space, 
with zero topological charge and whose total energy, total 
field momentum and z-component of the total magneti- 
zation, 

M z = J dx S z (x), 

are constants of motion. Writing the components of the 
spin operator as S x + iS v = Se l ^s'm8 and S z — Scosd, 
where S is the spin of the system, the general form of the 
dynamical soliton can be written as2^ 

e = 6(r-vt), (j) = wt + (p(r-vt). 

Here, v is the translational velocity of the soliton and 
w is the precessional frequency of the magnetization in 
the frame of reference moving with the soliton, i.e., an 
internal degree of freedom. 

A possible relation between the dynamical soliton and 
the bound states of n-magnons of the Heisenberg model 
was discussed by Schneider^ For the isotropic model, 
the dynamical soliton solutions of Tjon and Wright 32 
were semiclassically quantized via the Bohr-Sommerfeld- 
de Broglie condition. Following this procedure, the z- 



12 



component of the total magnetization assumes only in- 
teger values. Consequently, the precessional frequency 
w is also quantized^ For spin S = 1/2, it was found 
a correspondence between the n-magnons and the dy- 
namical soliton spectra. This result implies that the 
dynamical soliton of the one-dimensional ferromagnetic 
Heisenberg model can be considered as a bound state of 
n-magnons. This analysis was also applied to the easy- 
axis and anisotropic exchange Heisenberg models, but 
the correspondence between the two spectra was found 
only in the limit of large quantum numbers. 

Those characteristics are similar to the solutions of the 
Sine-Gordon models Two possible solutions of the clas- 
sical equations of motion are the topological soliton (a 
static localized solution) and the breather solution, which 
resembles an oscillating soliton-antisoliton pair (dynami- 
cal soliton) with topological charge equal to zero (a good 
review about topological and dynamical solitons is pre- 
sented in Ref. I3.1I) . After the quantization of these solu- 
tions, the soliton corresponds to a quasiparticle (fermion 
of the massive Thirring model) while the internal de- 
grees of freedom (oscillating modes) of the breather solu- 
tion correspond to bound states. Since the lowest energy 
bound state can be considered as an elementary boson of 
the theory, the internal degrees of freedom of the breather 
solution correspond to bound states of n-bosons. The 
number of the latter is determined by the coupling con- 
stant of the theory. 

It is well known that the low lying charged excita- 
tions of the two-dimensional electron gas at v = 1 is the 
quantum Hall skyrmion, which carries an unusual spin 
distribution. 16 As we will see in the next section, this 
excitation is described by a generalized nonlinear sigma 
model in terms of a unit vector field n(r) which is re- 
lated to the electronic spin [see Eq. (|62f> ] . The skyrmion 
is given by the topological soliton solution of the nonlin- 
ear sigma model with a finite size, which is determined by 
the competition between the Coulomb and Zeeman ener- 
gies. For v = 1, the topological charge of this solution is 
equal to the electrical charge. 

As our bosonization method for the 2DEG at v = 1 
gives us an interacting boson model to describe the in- 
teracting two-dimensional electron gas at v = 1, it seems 
reasonable to study the bound states of two-bosons in 
order to find out a possible relation between them and 
the spectrum of a bound skyrmion-antiskyrmion pair. 

Since the total and relative momenta of a boson pair 
are given by P = p + q and Q = (p — q)/2 respectively, 
the interacting term of the bosonic Hamiltonian (|50() can 
be written as 

«,„, » £ n^-Ktiv^f fc MP/2 + QA 

k,P,Q \ 1 J 

. / kA(P/2-Q) \ 
XSm i 2 ) 

X6 P/2+k-Q 6 P/2-k+Q 5 P/2+Q fc P/2-Q- (52) 



We can easily see that a state of two-bosons of the 
kind |$) = bVA\FM) is not an eigenstate of the total 
Hamiltonian (|50fl . Therefore, it is necessary to consider 
a linear combination of those states. The more general 
form of a state of two-bosons with total momentum P 
can be written as 

l$ P ) = $ p(q)^|p- q ^ P+q l™)- (53) 

q 

In this case, the total momentum of the pair is also a 
good quantum number for the same reasons as discussed 
at the end of Sec. Ill Al Remember that a state of two- 
bosons can be considered as a two electron-hole pairs 
whose total charge is zero. 

For a fixed value of the total momentum P, the energy 
of the state (|53() is given by the Schrodinger equation 

H|$ P ) =£p|$ P ). (54) 

We will consider the action on |<f>p) of the quadratic and 
quartic terms of the total Hamiltonian i|50|l separately. 
For Ho, we have 

Wo|*P> = E^^^'^P-q^P+q]!™) 
q 

+ J2^)b{ F _/ iF+ H \FM) 
q 

= X!*p(q)(wip_ q + wi P+q ) 

q V * ' 

Bp(q) 

x^P-ZiP+q™ 

+ E <f> p(q) & |p-q 6 lp + q^ M|FM> 
q 

= {Ep(q)+E FM )\<f>p). (55) 

Observe that £p(q) is the energy of two noninteracting 
bosons. On the other hand, for the Hint, after some 
algebra, it is possible to show that 

H mt |$p) = E^^^'^P-q^P + J 1 ^ 

q 

+E $ p(q) b iP-4p +q ^i FM ) 

q 

= 2 £ [/(k,P,q)$ P (q) 

k#0,q 

x^P-q+^W+q-^' (56) 

where 



Substituting expressions (|55|l and i|56|l in the Schrodinger 
equation 154|) and changing q — * q + k, we have 







5>p(q)(£ P (q) +E FM - E P )b{ b{ \FM) 



+ 2 C/(k,P,q)<i>p(q + k)6| p _ q 6t p+q |FM). 

k#0,q 

Changing the sum over momenta to an integral 

J 4tt 2 ' 

q 

we find the following eigenvalue problem 

(e-£p(q))$ P (q) = y"d 2 fc^ P (k-q,q)$p(k), (57) 

where e = Ep — Epm and the kernel of the integral equa- 
tion is given by 



*-p(k-q,q) 



Hk-q| 2 /2 

) _ | 

'tt |k-q| 

(k - q) A (P/2 + q) 

x sin 




(58) 



In the two expressions above, all momenta are measured 
in units of the inverse magnetic length, i.e., q — > q/Z. 

For the one-dimensional Heisenberg model, the analog 
eigenvalue problem can be solved analytically as the ker- 
nel of the integral equation is separable However, 
ifp(k — q, q) is not of the same kind and therefore our 
eigenvalue problem (|57|l will be solved numerically. 

The numerical solution of the above eigenvalue prob- 
lem can be determined using the quadrature technique. 29 
This method consists of replacing the integral over mo- 
mentum by a set of algebraic equations 

(e - E P (q i ))$p(q i ) « £ C^ P (qj - qi,q,)* P (qi)> 

(59) 

where Cj are the quadrature coefficients. The system of 
equations can be symmetrized multiplying then by \fCj, 

( e - J Bp(q i ))(^ /2 $ P (q i ))« 
^Cf 2 #p(qj - q i ,q i )C; /2 (C; /2 $ P (q j )). (60) 

After this discretization, for a fixed value of the total 
momentum P, we can calculate the eigenvalues of the 
equation H6()|l using usual matrix techniques. 

The choice of the points qj and of the values of 
the coefficients d are related to the parametrization 




FIG. 4: Dispersion relation of the state of two-bosons, in 
units of the Coulomb energy e 2 /(d) as a function of the total 
momentum P for g — 0. Solid line: lowest energy bound state; 
shaded area: continuum of the scattering states; dashed line: 
energy of the quasielectron-quasihole pair plus one spin wave 
[Eq. 161H 1 . See text for details. 



adopted. For one-dimensional problems, there are 
several parametrizations, for instance, the Gaussian 
quadrature^ 9 , which allows us to calculate the eigenvalues 
with good precision. However, for two-dimensional prob- 
lems, there are a fewer number of available parametriza- 
tions and therefore it is possible only to find a good esti- 
mate for the eigenvalues. 

In order to solve the eigenvalue problem (|6U|) we consid- 
ered a parametrization which is applied to calculate two- 
dimensional integrals over a circular region. In this case, 
it is necessary to introduce a cutoff for large momentum 
in order to define the integration region. All the parame- 
ters of this quadrature as well as a set of parametrizations 
for multiple integral calculations can be found in Rcf. EE 

Fig. 0] shows the dispersion relation of the states of 
two-bosons as a function of the total momentum P. Here, 
we assume that g = 0. The eigenvalue problem was 
solved using a 61 point quadrature and only bosons with 
momentum |kZ| < 2 were considered. The solid line is 
the lowest energy eigenvalue state while the shaded area 
is the continuum of scattering states. Once the low- 
est energy state of two-bosons is below the continuum 
of scattering states, we can say that this lowest energy 
eigenvalue (solid line) corresponds to a bound state of 
two-bosons. There are also other bound states above the 
one shown in Fig. but the analysis of those states is 
limited by the numerical method. 

As pointed out at the begining of this section, we want 
to check if there is a possible relation between the bound 
states of two-bosons and a bound skyrmion-antiskyrmion 
pair excitation. In this way, we should compare our re- 
sults with the ones derived from the model of Sondhi et 
alw&, namely, with the value of the energy of a nonin- 
teracting skyrmion-antiskyrmion pair, which can be cal- 
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culated from the expression derived by Sondhi et al. for 
the energy of the skyrmion [Eq. (7) of Ref 16]. However, 
this kind of comparison is not appropriate here as the 
z component of the total spin of the system is a good 
quantum number. Notice that the states of two-bosons 
has S z = 2 (the z component of the total spin in rela- 
tion to the quantum Hall ferromagnet) while a skyrmion- 
antiskyrmion pair described by the Sondhi's model has 
S* ^> 2. Remember that this model is suitable to de- 
scribe the skyrmion only in the limit of very small Zee- 
man energy (<? — > 0). In this case, the excitation is con- 
stituted by a large number of spin- flips. 

In a work previous tol&, Rezayi constructed a family 
of wave functions for the quasiparticles of the 2DEG at 
v = 1. Based on numerical calculations, it was shown 
that the energy (E\) of a state formed by the quantum 
Hall ferromagnet plus one spin down electron [see Fig. 
EI a)] is greater than the energy (E 2 ) of the state consti- 
tuted by one spin down electron plus a spin-wave exci- 
tation [Fig. Elk)]. In the thermodynamic limit, it was 
shown that Ei—E 2 = 0.054e 2 /eL This result implies that 
instead of a single spin down quasielectron, the quasipar- 
ticle of the 2DEG at v = 1 should be constituted by 
a quasielectron bound to n-spin waves. Based on that, 
Sondhi and coworkers suggested that the charged exci- 
tation of the 2DEG at v = 1 should be described by a 
charged spin texture. 

Notice that we can compare the spectrum of the bound 
states of two-bosons with Rezayi numerical results. 
Let us consider a state constituted by a quasielectron- 
quasihole pair very far apart and a spin wave with mo- 
mentum \lk\ -C 1 bound to either the quasielectron or the 
quasihole. The z component of the total spin of this state 
is S z — 2. Since the energy of a quasielectron-quasiholc 
pair very far apart corresponds to the limit \lk\ — > oo 
of the dispersion relation (|5T|l . the energy of the state 
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FIG. 5: Schematic representation of the quasiparticules con- 
sidered by Rezayi in Ref. l37l (a) the quantum Hall ferro- 
magnet plus one spin down electron and (b) the spin down 
electron plus a spin-wave excitation. 



describe above is simply 

Ee-h-™ « ( A?2 - 0.054)e 2 /eZ. (61) 

The dashed line on Fig. ^corresponds to E e -h- sw - 
We can see that our results are in good agreement with 
the previous ones of Rezayi's. More precisely, our calcula- 
tions indicated that, in the limit \IP\ — > oo, the dispersion 
relation of the bound states of two-bosons is asymptotic 
to E e ^h-sw In this scenario, we can understand the 
behavior of the dispersion relation of the bound states 
of two-bosons. As the total momentum \IP\ decreases, 
for instance, the quasielectron bound to the spin wave 
approaches the quasihole, increasing the interaction be- 
tween them and therefore lowering the energy of the sys- 
tem. Notice that this behavior is in good agreement with 
the solid line on Fig. Q] 

Therefore we can conclude that the bound states of 
two-bosons are appropriate to describe the skyrmion- 
antiskyrmion pair excitation, in the limit of large Zeeman 
energy, when the excitation is formed by a small number 
of spin- flips. 

We should mention that Cooper— studied the dynam- 
ical soliton solutions with zero topological charge of the 
nonlinear sigma model without the extra terms of the 
Sondhi's model. The calculated spectrum is qualitative 
similar to the one illustrated in Fig. 0] For small mo- 
mentum, the excitations correspond to free spin waves 
and, as the momentum increases, the dispersion relation 
continuously approaches the energy value of a noninter- 
acting soliton-antisoliton pair. However, this description 
is valid only in the limit of large number of spin waves, 
which is very far from the region of our analysis. 

A final word about the Hilbert space. Notice that in 
the above analysis we consider the bosonic Hilbert space 
constituted by the vacuum state [Eq. l(TU|l ] and the states 
of two-bosons [Ed. 1(53)1]. As discussed in Sec lIIEl the 
number of the states of two-boson is greater than the 
number of fermionic states with two spin-flips over the 
quantum Hall ferromagnet. However, in order to solve 
the eigenvalue problem l|60|) , a cutoff for large boson mo- 
mentum was introduced which restricted the bosonic sub- 
space and therefore we believe that the solution of Eq. 
(|60|1 does not involve unphysical states. 



IV. BOSONIZATION AND COHERENT STATES 

In this section, we will consider the semiclassical limit 
of the interacting bosonic Hamiltonian (|50|l . We will 
show that, starting from l]5Up. it is possible to recover 
the energy functional of the quantum Hall skyrmion. 

As mentioned in Sec. IIII CI Sondhi et alM suggested 
that the quantum Hall skyrmion can be described by a 
generalized nonlinear sigma model in terms of an unit 
vector field n(r) which is related to the electronic spin. 
The effective Lagrangean density of the model is given 
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by 



£ e // = ^p A(n) ■ d t n- ips(Vn) 2 + \g*p p, B n- B 



2e ./ |r - r' 



(62) 



Here, ps is the spin stiffness (see Appendix lA")l . A(n) is 
the vector potential of an unit monopole (e ahc d a Ab = n c ), 
po = 1/ (27rZ 2 ) is the average electronic density, and q(r) is 
the topological charge density or skyrmion density which 
is given by 



q(r) = —e aP n ■ (d a n x dpn), 

07T 



(63) 



with a,b,c — x, y, z, a, [3 — x, y, and e a @ is the antisym- 
metric tensor. 

On the other hand, Moon and co-workers suggested 
an alternative approach to study the quantum Hall 
skyrmionii& In this case, the charged excitation is de- 
scribed by the state 



|n(r)> = e~ i0 \FM), 



(64) 



where the operator O is a nonuniform spin rotation which 
reorients the local spin from the direction z to n(r), 

O = J d 2 r J2(r) ■ S(r) 

= J d 2 q [n-(q)sl q + n+(q)si q ] 



(65) 



Here, S(r) is the spin operator, n(r) is a unit vector, 
and fi(r) = z x n(r) defines the rotation angle. Assum- 
ing that O(r) corresponds to small tilts away from the z 
direction, S! CT (q) vanishes when \lq\ 3> 1. 

In this long wavelength approximation, it was shown 
that the average value of the electron density operator in 
the state (|6^)l is equal to the topological charge density 
of the vector field n(r) [Eq. (J53J]. Moreover, after pro- 
jecting the Coulomb potential in the lowest Landau level 
subspace, its average value in the state (|64|l is equal to 
functional energy derived from the Lagrangean density 

In Sec. IIII CI we showed a possible relation between 
the skyrmion-antiskyrmion pair excitation and the bound 
states of two bosons. Therefore if we consider a semiclas- 
sical limit of the interacting bosonic Hamiltonian (|50|l in 
the same way as it was done in Ref. 18: we can check if 
it is possible to recover the results of Sondhi et ali&> 

Let us assume that equation (|64|l is a good approxima- 
tion to describe the skyrmion. Substituting the expres- 
sions l|31[) and (|29[l in (|65|) and approximating S+ only 
by the linear term in the bosonic annihilation operator, 
we can write down the state (|64[l as a function of the 
bosonic operator b as 



\sk) 



e-"e- i0 



\FM), 



(66) 



where the operator O is redefined as 



(67) 



and the constant M = 2 ^^ 2 J d 2 q J7+Sl_ q . The value 
of the constant (3 will be determined later. Observe that 
the state (|66fl is a coherent state of the bosons b. 

Changing the sum over momenta into an integral in the 
expression of the bosonic representation of the electron 
density operator /5k [Eq. I|27|l] the average value of this 
operator in the state (|66[) is given by 



(sk\p\\sk) 



2 (2tt) 5 

x J d 2 q sin(k A q/2)n+fi q+k . (68) 

In the above expression, the momenta are measured in 
units of the inverse magnetic length /. In the long wave- 
length approximation (remember that ri CT (q) is different 
from zero only when \lq\ <C 1), we have 

e- fc2/4 sin(k A q/2) w k A q/2 = z ■ (k x q)/2, 
and therefore Eq. Q68f) can be written as 

(sk\ Pk \sk) = A^((q + k)O q+k ) x (qO+) 



2 5 7T 3 



e afi I d 2 r e- lk r z • (Vn Q x Vr/) 



(69) 



In the second step, we use the fact that fi(r) = z x n(r), 
hence Q x = —n y , Q y = n x and fl z = 0. From equation 
(|69|l . the Fourier transform of /5k is given by 



p(r) = (sk\p r \sk) 
1 . 



4tt 2 8tt 



e Q/3 i • (Vn a x V? 



(70) 



As we have assumed that fi(r) corresponds to a small 
rotation angle of the local spin, Vn z w and n z ps 1. 
Therefore, we can write 



e a pz ■ (Vn Q x Vn 13 ) w e a pn z z ■ 



(Vn Q x Vr 



w e ap n ■ (9 Q n x dpn) , 

as it was done in Rcf. Moreover, if we choose the 
value of the constant = Att 2 /N 2 , Eq. JJJJ is in agree- 
ment with the definition of the topological charge density 

Following the same approximations, we will calculate 
the average value of the energy of the state \sk) for the 
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interacting bosonic Hamiltonian l|50|l . The average value 
of the quadratic term of the Hamiltonian (|50|l is given by 



(sk\H \sk) 



4(2 



4)3 / A- q O± q O q . (71) 



Considering the long wavelength limit of the dispersion 
relation l|51l) . i.e, 



1 



w q ~ 9 + ^b\ 



Eq. I|71|l can be written as 

(sk\H \sk) « (H z } + {Hg), (72) 
where the Zeeman term is given by 

{7iz) = 4^ 9 J d2qn+ -^ 

and the gradient term by 

(H S ) = f 4W/^^ 2 ^- 

Here, the constant e# is defined in Appendix 1X1 

Rescaling the momenta by Z" 1 and calculating the 
Fourier transform, we can show that the Zeeman term 
can be written as 

(sk\n z \sk) = \g*»BB^ J ' d 2 r(jff + { n yf 

1 1 f 

« ~2 5 * MB 27 / U 

+ \9>bB^-N^. (73) 

Z Z7T 

In the second step above, we use the identity 

|n - z\ 2 = {n x f + {n y f + (n z - l) 2 = 2 - 2n ■ z 



(1201, we have 

{'Hint) = 



2ANj 



-|ifc| 2 /2 



sin (k A p/2) 



k,p,q 



sin(kAq/2)r! k+p r! q _ k Ot p fi: q 
' f d 2 kV(k/l) 



2Z 2 (8tt 2 ) 2 

< | d 2 pz-((p + k)r! p+k )x(pr!+ 

< | A^((q-k)fi-_ k )x(qO+ 
- J d 2 rd 2 r' V{r-r')q{r)q{r'), 



(75) 



where V(r — r') = e 2 /(e|r — r'|) is the Coulomb potential, 
q(r) is the topological charge density as defined in [Eq. 
Q70[l] and the vector r is measured in units of /. 

From Eqs. (J73J, l(7^|) . and J7SJ), we can conclude that 
the average value of the energy of the state \sk) is given 
by 



(ak\H\8k) = -p% 



d 2 r [Vn(r)f 



Z Z7T 



2 9 ""fr 



d 2 r n • B 



+ 



2 .2 / g( r )g( r/ ) 

, I u rd r — — . 

2d. r-r' 



(76) 



Notice that Eq. (|76|) is equal to the energy functional 
derived from the Lagrangean density Q62JI for a static 
configuration of the vector field n(r). 

It is important to mention that the choice of the con- 
stant f3, based on the fact that Eq. H7()|) should be equal 
to the topological charge density l|63|l related to the vec- 
tor field n(r), gave us the correct values of the coefficients 
of the Zeeman (|73|) . gradient l|74() . and interacting l|75() 
terms. 

In addition to the approach discussed in Rcf. [l8|, a 
different method to derive an effective field theory for 
the quantum Hall skyrmion is presented i n ^2ii2i2i 



and the fact that, within our approximation, \n z — 1| <^ 1. 
On the other hand, the gradient term assumes the form 

( S k\H \sk) = ip s J d 2 r [(V^) 2 + (V^) 2 ] 

« \ps J d 2 r [Vn(r)] 2 , (74) 

where ps is the spin stiffness as defined in Eq. {Jj2J|. 
Finally, for the interacting term of the Hamiltonian 



V. SUMMARY 

We developed a bosonization approach for the 2DEG 
at v = 1 using the fact that, at some level of approxi- 
mation, the elementary neutral excitations of the system 
can be treated as bosons. The Hamiltonian of the 2DEG 
at v = 1, the electron density, and spin density operators 
were bosonized. We showed that the bosonic representa- 
tion of the spin density operators is analogous to the one 
considered by Dyson to study the ferromagnetic Hcisen- 
berg model. Furthermore, we showed that the developed 
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bosonization method is closely related to the LLL pro- 
jection formalism developed by Girvin and Jach. 

The method was applied to study the interacting two- 
dimensional electron gas at v = 1. The Hamiltonian 
of the fermionic system was recast in an interacting two- 
dimensional boson model. We showed that the dispersion 
relation of the bosons is equal to the previous diagram- 
matic calculations of Kallin and Halperin. Within our 
bosonization approach, we can go beyond the latter re- 
sults as we also found an interaction between the bosons. 

Finally, we showed that the derived interacting bosonic 
model can describe the quasiparticle-quasihole pair exci- 
tation of the 2DEG at v = 1. On one hand, we showed 
that the interaction between the bosons accounts for the 
formation of bound states of two bosons. Our results 
agree with the previously developed numerical approach 
of Rezayi's, who studied the quasiparticles of the sys- 
tem. On the other hand, we showed that the semiclassi- 
cal limit of the interacting bosonic Hamiltonian recovers 
the energy functional derived from the model suggested 
by Sonhdi et al. to describe the quantum Hall skyrmion. 
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APPENDIX A: ENERGY AND LENGTH SCALES 
FOR THE 2DEG 



In the table below, we show the cyclotron Hw Cl Zeeman 
g and Coulomb e c energies and the value of the constants 
cb and p s , in Kelvin, as a function of the magnetic field 
B. The magnetic length I = ^hc/(eB) = 256/VB is 
measured in angstroms and the magnetic field B in Tesla. 



Energy scales 


(K) 


hw c 


heB/(m*c) 


18.785 


9 


g*H B B 


0.33S 


€c 


e 2 /(d) 


50.40/B 


es 




63.16%/^ 


Ps 


e c /(16V2^) 


1.25VB 



The electron effective mass in the GaAs quantum well 
is to* = 0.0777i e , where m e is the electron mass and the 
dielectric constant of the semiconductor is e w 13. 



APPENDIX B: CHARGED PARTICLE IN A 
PERPENDICULAR MAGNETIC FIELD 

Let's consider an electron moving in the x — y plane 
under the action of a constant perpendicular magnetic 
field B = Bqz. The Hamiltonian of the system is given 
by 



H = 



-^(p + -A(r)) 2 

2m V c / 



(Bl) 



where p is the momentum canonically conjugated to r 
and A is the vector potential. In the symmetric gauge, 



A(r) 



1 



r x B 



1 



B (yx-xy). 



Classically, the electron moves in a circular orbit with 
angular frequency w c — eB/mc (cyclotron frequency). 
In this case, the modulus of the particle velocity 



v = — fp + - A(r) 

m V c 



and the position of the center of the cyclotron orbit (guid- 
ing center) 



Ro = r- 



z x v 



are constants of motion. 

Defining the complex variables V = v x + iv v , P = 
Px + ipy, Z = x + iy, and Zq — Rq x + iRo v , in the 
symmetric gauge, the constants of motion can be written 
as 

V = -P+\iw c Z, 

777 2 

Z Q = Z+—V. 

w c 

Now, if we apply the canonical quantization rule for the 
canonical conjugate variables r and p, the commutation 
relations between V and Zq are given by 



[V,V*] = - 



2hw r 



-2l 2 w 2 c , 



2i 2 



[Zo,Z Q \ — 

777 W c 

[VKZ Q ] - [V,Z o ]=0, 



where I is the magnetic length. Introducing two indepen- 
dent ladder operators d and g, such that [d, d<] — [g, g'] — 
1 and [d, g] = [d, g'] = 0, we can write 



V = -iV2lw c d\ 



(B2) 



Zq = V21 



fl- 



it is easy to prove that the operators V and Z defined 
as in Eq. I|B2(1 satisfy the above commutation relations. 
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Therefore the Hamiltonian (IB II) can be written as 



Ha = hw c (Sd + -), 



(B3) 



whose energy eigenvalues (Landau levels) are given by 



E ntm = hw c (n + -) 



(B4) 



and the energy eigenvectors by 



(r|0,0) = 
(r\n m) = 



1 _, 



10 0), 



1 



« G m+n , n ( t ). 



(B5) 



Here the function G m + n . n (x) is dehncd in Appendix IO 
Semiclassically, the state \n m) can be seen as an elec- 
tron in a cyclotron orbit with radius equal to l\/2n + 1 



and the center located at a distance ly/2m + 1 from the 
origin of the coordinate system. 

A detailed analysis of this problem is presented in Ref. 
Our formalism is similar to the one presented in this 
reference with the replacements g — > w and <v — ► — ia. 



APPENDIX C: THE G m . m ,(lq) FUNCTION 
PROPERTIES 



We want to calculate the matrix element of the opera- 
tor e~ 4qr in the Landau level basis. Writing q — q x + iq y 
and r = x + iy, we can expand the latter in terms of the 
ladder operators d and g defined in Appendix iBl [see Eq. 

Ell, 



r = Z = Z V = V2l(g - d f ), (CI) 



Therefore the matrix element becomes 



(nm\e lq ' r |n'm') = (nm\ exp (— i(qr* + q*r)/2) |n'm'} 



{nm \ exp 



-il ((qg + q*g^) - (q*d + qttf)) /V^] \n' to') 



(C2) 



Since the ladder operators d and g are related only to the Landau levels and the guiding centers respectively, we can 
use the properties of these operators to write the above matrix element as a product 



(nm\e-^ r \n' m') ^ eM~H 2 /2)G m , m ,(lq)G n ,A-lq*), 
where the functions G m ^ m i(lq) and G n>n i(—lq*) are defined as 

G m ,m'(lq) = (m\ exp(-j^5 t /v / 2) exp(-ilq*g/y/2)\m'), 



(C3) 
(C4) 



G n , n -{-lq*) = (n\exp(ilq*d/V2)exp(ilqdyV2)\n'). 

Now, if we take n = n' = in Eq. I|C3() . we have the matrix element of the operator e~ iqr in the lowest Landau level 
basis 



He-^lm') = exp(-\lq\ 2 /2)G m , m .(lq). 



(C5) 

— iq-r 



In the LLL projection formalism, the above expression corresponds to the matrix element of projected operator e 
[compare Eq. with Eq. (25.1.11) of Ref. US$). 

Using the properties of the ladder operators, it is possible to show that the function G m ^ m '(lq) can be written as a 
linear combination of the generalized Laguerre polynomials L™7 m (|/<7| 2 /2), i.e., 



m'\ V \/2 

(m — m)\ r ■=■ 

V ml \ y/2 



V 2 



j m— 



(C6) 



From expressions (|C4|) and l|C6(l we can prove the following properties of the function G m>m j (Iq) 
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(i) relations between the function and its complex conjugate: 

G m ,m'{lq) = G mrn i(—lq ) = G^t >m (—lq) = G m ^ m (lq*) 

G m , m >(ilq) = G* m<m ,(ilq*) = G* m , tm (-ilq*) = (-i) m ~ m ' G m> , m (lq*). (C7) 

(ii) The Fourier transform of the product of two functions: 

e^G m ,m'(lq)G n , n ,{-lq*) = J d 2 r e~ l * r (n', m'|r) (r\n, m) 

= i J d2r e " iq re ^ G »+™.« (y) G «',™'+«' (^) • (C8) 

(in) The sum of the product of two functions: as the Landau level basis \n, m) is a complete basis, we have 

G mJ {lq)G Um > (Ik) = 53(m| eM-Hqtf /V2) cxp(-ilq*b/V2)\l) (l\ exp(-ilktf /V2) exp(-ilk*b/V2)\m') 
i i 

= exp (^-£^) G m , m ,(lq + Ik). (C9) 



(iv) Orthogonality relation: using the orthogonality relations of the generalized Laguerre polynomials, we can show 
that 

d 2 k e -|^| 2 /2 Gm m , (-lk*)G n , n . (Ik) = ^6 m , n S m >,n> , (CIO) 

and changing the integral over momenta by a sum, 

J2 e- m2/2 G m , m ,(-lk*)G ntn ,(lk) = AVW-W- (CH) 

k 

(v) The trace: 

P l'<?l 2 /2 r _ k i2 j r 

Y,G m , m (lq) = d 2 re~^e^G Q , m ( — )G m ,o(j) 

m m 

e |i<?l 2 /2 r _, r ,2 | H 2 

= / e-'^ r e^ e w G ,o(0) - A^(q). (C12) 



APPENDIX D: THE COMMUTATOR [S+,5~,] 

If we consider the expressions of the spin operators S£ and S~ in terms of the fermionic annihilation and creation 
operators [Eq. l(T^|l and iJEJ], we have 

K,S-,] = e -M 2 /*-W\ 2 /i £ G^A^G^Wttci^^c^c,^} 

= e -\l q ?l2-\l q '\-l2 I J- G m ,„Gg)G nX (/g')4T C «'T- E G »,m(««')G ! m,m'(/g)4 ; C roU 
\ m,n,n' m,m' ,n 

\ m,n m.n / 

= e^'^e-'^'' 2 / 2 ^^^ (D1) 
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Now, if we compare the above result with the expressions of the electron density operators /5 CT (q) [Eq. (TTJ], we can 
conclude that 

K, S-,] =e^'*/ 2 p T (q + qO-e^*/^ ;(q + q '). (D2 ) 

I 
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